NN SCATTERING: CHIRAL PREDICTIONS FOR 
ASYMPTOTIC OBSERVABLES 

J-L. Ballot 1 * 

Division de Physique Theorique, Institut de Physique Nucleaire 2 ^ 
F-91406, Orsay CEDEX, France 



M.R. Robilotta 3 * 
FINPE, Instituto de Fisica, Universidade de Sao Paulo 
CP. 66318, 05389-970, Sao Paulo, SP, Brasil 

C.A. da Rocha 4§ 
Department of Physics, University of Washington 
Box 351560, Seattle, Washington 98195-1560, U.S.A. 
(March, 1997) 



Abstract 

We assume that the nuclear potential for distances larger than 2.5 fm is given 
just by the exchanges of one and two pions and, for the latter, we adopt a 
model based on chiral symmetry and subthreshold pion-nucleon amplitudes, 
which contains no free parameters. The predictions produced by this model 
for nucleon-nucleon observables are calculated and shown to agree well with 
both experiment and those due to phenomenological potentials. 
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I. INTRODUCTION 



In the last two decades several semiphenomenological potentials were proposed [l]-[5[ 
which can describe rather well NN scattering data below the pion production threshold. 
All these models have a common feature, namely that the interaction at large distances 
is ascribed to the one-pion exchange potential (OPEP). On the other hand, they differ 
significantly at intermediate and short distances, a fact that becomes evident when one 
inspects the profile functions produced for the various components of the force. In each case, 
the reproduction of experimental data is achieved by means of a different balance between 
effects due to long and short distances. The fact that these potentials are successful means 
that somehow they are able to incorporate the relevant average dynamics. 

At present, none of the existing semiphenomenological potentials include explicitly the 
dynamics associated with chiral symmetry, which constitutes the main conceptual framework 
for the study of strong interactions at energies which are small compared to the QCD 
scale. In this regime, non-perturbative effects are dominant and one is not able to do 
calculations using QCD directly. The usual strategy for overcoming this problem consists 
in working with an effective theory, constructed in such a way as to include as much as 
possible the main features of the basic theory. The masses of the quarks u and d are very 
small and hence their interactions with gluons are approximately invariant under the group 
SU(2) x SU(2). Therefore, one requires the effective theory at the hadron level posses the 
same basic symmetry, broken by just the pion mass. 

In the last five years several authors have tackled the problem of NN interactions in the 
light of chiral symmetry and, so far, only processes associated with the two-pion exchange 
potential (TPEP) were systematically studied p-^. Chiral symmetry is very relevant 
to this component of the force because it controls the behaviour of the intermediate 7rN 
amplitude, that is the main building block of the interaction. The first works of this series 
were restricted to systems containing just pions and nucleons and considered basically the 
first five processes given in Fig.l (A) |6|- [l0l , [T2| . These processes constitute an autonomous 
chiral family, incorporate correctly the well known cancellations of the intermediate 7rN 
amplitude [16|jl7]] but correspond to an intermediate amplitude which is too simple for 
reproducing 7rN experimental data fl8fl . 



The extension of this minimal model so as to include other degrees of freedom was con- 
sidered by Ordonez, Ray and van Kolck (|TT|,|T^] and by ourselves [ |H5| , |H)| | . In the former case, 
a very general effective Lagrangian was used, which included explicitly the interactions of 
pions, nucleons and deltas and contained some free parameters representing other interac- 
tions. In principle, these parameters could be obtained from other physical processes, but 
these authors chose to adjust them to NN scattering data. Using a non-relativistic cut-off of 
the order of the rho mass, they could achieve a qualitative description of all NN observables. 

In our approach, the intermediate 7rN amplitude involved the interactions of pions and 
nucleons, supplemented by empirical information in the form of th Hohler, Jacob and Strauss 
(HJS) coefficients |19|]. In general, the physical amplitude for the process T[ a {k)N{p) — > 
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Tr^(k')N{p') may be described by two independent variables, v = ^(p + p') ■ (k + k') and 
t = (k — k') 2 . In order to obtain the HJS coefficients, one subtracts the nucleon pole 
from the empirical 7rN amplitude and uses dispersion relations to continue analytically the 
remainder (R) to an unphysical region around the point v — 0, t — 0. The HJS coefficients 
are then obtained by expanding this remainder in a power series in v and t. The use of 
the subthreshold coefficients is particularly suited to the calculation of the TPEP at large 
distances, since this part of the potential is determined by the intermediate 7rN amplitude 
in the very neighbourhood of the point v = 0, t = 0. Thus, in this aspect, our approach is 
very different from that of Ordonez, Ray, and van Kolck. 

In our calculation, the TPEP was derived from the ten diagrams depicted in Fig.l, 
representing both interactions involving only nucleons (A) and other degrees of freedom 
(B). The main features of the asymptotic TPEP were extensively discussed in jnj and here 
we just recall some of the conclusions of that paper. One of them is that the scalar-isoscalar 
component of the TPEP at large distances is attractive and therefore qualitatively coherent 
with a well known feature of the nuclear force. As far as dynamics is concerned, we have 
shown that chiral symmetry is responsible for large cancellations in the pure nucleon sector 
(Fig.l, A) (l^ll^l and eventually the contributions from this sector turn out to be much 
smaller than those arising from other degrees of freedom. The main contribution to the 
intermediate attraction is due to processes containing nucleons in one leg and the remainder 
degrees of freedom in the other. 

The fact that our calculation of the TPEP did not contain free parameters means that 
it yields predictions for NN observables, whose study is the main goal of the present work. 
We assume that, for distances larger than 2.5 fm, the NN interaction is given by just the 
OPEP and the TPEP calculated in ref. fl5| and try to determine the values of the angular 
momentum and the energy regions for which observables can be ascribed to these components 
only. Our presentation is divided as follows: in Sec. II we discuss our method of work and 
in Sec. Ill we give our results and conclusions. 



II. DYNAMICS 

In general, it is not easy to isolate unambiguously the observables associated with a 
particular region of a given potential. Nevertheless, in many cases, the centrifugal barrier 
can suppress a significant part of the short range interaction and one is left only with 
contributions from the tail of the force. For instance, in a study of the influence of the 
OPEP over NN observables, we have obtained, for most waves with i > 2, purely pionic 
phase shifts and mixing parameters, which did not depend on the short range features of 



the interaction ||20|| . In the case of the OPEP, this kind of result is possible because the 
potential is not too strong. 

In the present problem, even if one is willing to consider the influence of the TPEP only 



for separations larger than 2.5 fm, where the results of Ref. [15] are mathematically reliable, 



one needs to use expressions which are valid for all distances. As it is important have close 
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control of the regions of the potential that contribute to the observables, we employ the so 
called variable phase method. It is fully equivalent to the Schrodinger equation and provides 
a clear spatial picture of the way phase shifts and mixing parameters are structured. For 
the sake of completeness, we summarize here the main equations used in our calculation. 
In the case of uncoupled channels, the wave function uj(r) with angular momentum J is 
written as [plfl 

Uj(r) = cj(k, r) jj(kr) - sj(k, r) fij(kr) , (1) 



where jj and fij are the usual Bessel and Neumann functions multiplied by kr. The functions 
cj and sj, for a potential Vj(r), are given by 

cj(k, r) = 1 - - j o dp Vj(p) jj(kp) uj{p) , (2) 



sj(k,r) = — - / dp Vj(p) hj(kp) uj(p) . (3) 
k Jo 

The variable phase Dj(k,r) is defined as 

Dj = tan~\ S -±) (4) 

and, by construction, it vanishes at the origin and yields the observable phase shift 5j when 
r tends to infinity. Differentiating (2) and (3), using (1), and manipulating the result, one 
obtains the differential equation 

D'j = -j Vj P 2 j(D 3 ) , (5) 

where D'j = and the structure function Pj is given by 

Pj = jj cos(Dj) - hj sin(Dj) . (6) 



In the case of coupled channels, one has two phases, Dj m and Dj p , with m = L — 1 and 
p — L + 1, and a mixing parameter Ej, that depend on r and become the observables 5j m , 
5j p and e,j when r tends to infinity. Denoting the diagonal and tensor components of the 



potential by Wjl and Tj, one has the following coupled differential equations |22 

[w. 



D' 



in 



Jm 



Jm 



cos\Ej) Pi - sm\Ej) Q 



k cos{2Ej 
W Jp sin\Ej) cos 2 (Ej) (P* - Q 2 p ) 

2 Tj sin(Ej) cos(Ej) cos 2 (Ej) P m Q p - sin 2 (Ej) P p Q m } 



(7) 
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E' Jm = - — {Tj [cos 2 {Ej) P m P p + sm 2 (Ej) Q m Q % 

- W Jm sin(Ej) cos(Ej) P m Q m - W Jp sin(Ej) cos(Ej) P p Q p } . (8) 

In these expressions the structure functions Pl and Ql are defined as 

Pl = 3l cos(Dj L ) - h L sin{D JL ) , (9) 
Ql = 3l sin(D JL ) + n L cos(D JL ) . (10) 



The equation for Dj p is obtained by exchanging the labels m and p in (7). 

As far as the interaction is concerned, we consider just the OPEP (V^) and the TPEP, 
which are assumed to represent the full potential for distances greater than 2.5 fm. As 
pointed out in the introduction, the TPEP is determined by two kinds of contributions, one 
generated in the pure pion-nucleon sector (Vn) and the other associated with the remainder 
degrees of freedom (Vr). The direct inspection of these potentials indicates that the former 
is comparable to the OPEP, whereas the latter is rather strong for r > 2.0 fm. In order to 
calculate the observables, one needs to regularize the various potentials at short distances. 
In the case of OPEP, the regularization is achieved by cutting it at a radius r n and replacing 
the inner part by the constant value Kr(fv)- As Vn is comparable to the OPEP, we adopt 
the same regularization procedure for it, with a radius r^. 

The regularization of Vr is more problematic. For distances around 1.0 fm, the value of 
the cental component of the potential is about —25 GeV. On the other hand, in Ref. [15] we 
have argued that the asymptotic TPEP is mathematically reliable only for distances larger 
than 2.5 fm, indicating that the odd behaviour of the TPEP at short distances is unphysical 
and associated with the use of equations outside their domain of validity. Inspecting the 
equations used in that work, it is easy to relate this behaviour to the HJS coefficients 
involving high powers of v and t in the intermediate 7rN amplitude. On the other hand, 
restricting ourselves to the first two leading contributions, due to the coefficients of the terms 
v°t° and z/°t 1 , we keep most of the asymptotic potential and get values around —150 MeV 
in the neighbourhood of 1.0 fm, which are still large, but much more reasonable. Therefore 
we base the present study on this leading potential, which is regularized by means of a step 
function 9{tr). 

In this calculation one can rely only on those results which are independent of the 
radii used in the regularization procedure. In order to control this independence, we adopt 
r^v = r n , vary in the interval 0.8 — 1.0 fm and discard the cases where the contribution of 
Kr + Vn to the observables varies more than 1.0%. Concerning Vr, the preceding discussion 
suggests that one should be interest in effects due to the region r > 2.5 fm and hence we 
consider values of tr between 1.5 fm and 2.5 fm and study the effect of this variation over 
the observables. This produces an indication of both the stability of the results and the 
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importance of the inner part of the potential. When the variation of the results is less than 
5%, we take them as predictions of the potential. For larger deviations, we consider them 
as estimates. 



III. RESULTS AND CONCLUSIONS 

We have calculated the predictions for NN observables produced by a chiral potential 
15] involving only the exchanges of one and two pions, assumed to represent the full inter- 
action for distances larger than 2.5 fm. Since we are interested only in the cases where the 
centrifugal barrier cuts naturally the inner parts of the force, we used the variable phase 
method to control this aspect of the problem. One may acquire a feeling for this method 
by looking at Fig. 2, where we display the variable phases for the wave 1 G±, divided by 
their asymptotic values, for three different energies. It is possible to see that, as expected, 
higher energies probe more the interior of the potential. It is also interesting to note that 
the variable phase method allows one to make quantitative statements such as, for instance, 
the radius for which the phase attains a given percentage of its final value. 

In the case of uncoupled waves, the interplay between the centrifugal barrier and the 
regularization of the background is rather intuitive, but this is not the case of coupled 
systems. In order to clarify this point, we show in Fig. 3 the variable phases for the 
3 D 3 — 3 G 3 system for E = 200 MeV, due just to the background and regularized at either 
at 0.7 fm or 1.0 fm. One sees that the 3 D 3 curves depend strongly on the cutoff, whereas 
those describing e 3 and 3 G 3 are very stable and cannot be distinguished with the naked eye. 
The wave 3 G 3 is negligible for distances smaller than 2.0 fm due to the centrifugal barrier, 
meaning that the system is effectively uncoupled up to that distance. The construction of 
3 D 3 extends up to 3.5 fm, where £3 is maximum, but the two other observables become 
asymptotic much later, around 5.0 fm. This indicates that, even for coupled waves, the 
various components get their contributions from different regions in space. 

We found out that the observables associated with the following waves depend more than 
1% on the cutoff used for the V n + Vn background and hence are not suited for our study: 
%, 3 P , 3 P 2 , 3 S 1 ,e 1 3 D 1 and 3 D 3 . 

In Figs. 4-7, we display our results for the phase shifts as a function of the laboratory 
energy. They include predictions from the OPEP cut at 1.0 fm, the sum of the OPEP and 
Vn cut at 1.0 fm and the total results produced by cutting the contribution of Vr at 1.5 
fm (x(1.5)) and 2.5 fm (x(2.5)), and the experimental values are taken from SAID VZ40 
solution [^|. For comparative purposes, we also include the predictions of the Argonne [§] 
phenomenological potential. 

The dominant part of the chiral TPEP is associated with the exchange of a scalar- 
isoscalar system and hence its significance depends strongly on the NN channel considered. 
Therefore in the sequence we comment the main features of our results for the various 
subspaces of spin and isospin. 

(T=1,S=0); Fig. 4: There are no results for the wave 1 Sq, since it cannot be understood 
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as a TPEP superimposed to a pionic background. For the waves X D 2 and 1 G± one has 
predictions for energies up to 50 Mev and 300 MeV respectively. As expected, the TPEP 
increases the attraction due to the OPEP and our results are very close to the experimental 
values for the potential cut at 2.5 fm. 

(T=1,S=1); Figs. 5a,b,c,d: The waves 3 Pq and 3 P 2 depend on the pion cutoff and 
were discarded. For the uncoupled waves 3 Pi, 3 P3, and 3 H 5 (Fig. 5a) we obtain predictions 
which extend up to 300 MeV for the last two of them. Results for the waves 3 Pi and 3 F% 
are compatible with experiment, but this does not happen for the 3 H 5 wave. In the case 
of the coupled waves, that with lowest orbital angular momentum tend, as expected, to be 
much more influenced by the cutoff used for the OPEP than the other ones. We obtain 
predictions in all cases, but the mixing parameters are heavily dominated by the OPEP and 
yield very little information about the TPEP. For the waves 3 F 2 (Fig. 5b) and 3 F 4 (Fig. 
5c), the differences with experimental values are small, whereas for the waves 3 H 4 (Fig. 5c) 
and 3 Hq (Fig. 5d) they are important. 

(T=0,S=0); Fig. 6: Our calculation yields predictions for all the waves in this sector, 
namely l P\, 1 i 7 3, and l H^ n generally quite close to the pure OPEP ones, reflecting the fact 
that our central potential in this channel is small. Results are also close to experiment. 

(T=0,S=1); Figs. 7a,b,c: The observables 3 Si, ei, 3 D 1 and 3 D 3 depend on the pion 
cutoff and are not considered. In all cases, coupled and uncoupled, our results are dominated 
by the OPEP and close to experiment. 

In order to assess the general trends of the various observables presented in Figs. 4-7, it 
is useful to recall that the relative strengthes of the central and tensor OPEP in the channels 
(T,S) = (1,0), (1,1), (0,0), and (0,1) are respectively 1 : § : 3 : 1 and : § : : 1. This 
means that one pion exchange is more important in the channels with T = and hence the 
good agreement between predictions and experimental results noted in Figs. 6 and 7 may 
be ascribed to OPEP physics. 

As far as the channels with T = 1 are concerned (Figs. 4 and 5), the tensor interaction 
makes the OPEP to be more important for triplet waves, in spite of its weaker central 
component. Therefore the role of the TPEP is more evident in the waves l D 2 and 1 G*4 
(Fig. 4), where chiral predictions agree well with experiment. In the case of triplet waves 
(Figs. 5a,b,c,d), one also finds that the chiral potential is able to reproduce experimental 
data when I < 5, but this does not happen for H waves. The behavior of these waves are 
peculiar, since they have a high orbital angular momentum and hence should be close of being 
OPEP dominated. Indeed, our results show that predictions from the chiral potential for H 
waves are not far from the OPEP and also depend little on the cutting radius. Part of the 
discrepancies observed may be associated with the fact that we have used g 2 /47r = 14.28 for 
the 7tN coupling constant |I{J whereas the SAID analysis is based on the value g 2 /4:7i = 13.7. 
It is also worth pointing out that there is a 10% difference between the experimental pn and 
pp solutions and the discrepancies would be reduced if the latter were used. However, even 
if these factors were considered, the experimental data would still seem to suggest that the 
TPEP is repulsive for the waves 3 Hi and 3 H$, something which is rather difficult to explain 
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theoretically. 

A general conclusion that can be drawn from the present study concerns the details of 
the TPEP. As discussed in the introduction and represented in Fig. 1, it consists in a sum 
of terms, arising from both the pure pion nucleon sector and from interactions involving 
other degrees of freedom. Our results show that the former contributions are very small, 
indicating that the numerical significance of the TPEP is essentially due to the interplay 
between nucleon and other degrees of freedom. 

In this work we tested the chiral TPEP derived in Ref. [jHJ, which is based on subthresh- 
old 7rN data and contains no free parameters. Our results have shown that it is rather 
consistent with experiment^ 
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Note added in revision: In a recent work, the predictions from a similar chiral potential were 
presented |M| , which agree qualitatively with those produced here. 
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FIGURES 



(A) 




FIG. 1. Contributions to the two-pion exchange nucleon-nucleon potential from the pure 
pion-nucleon sector (A) and for processes involving other degrees of freedom (B). The crosses 
in the nucleon lines of the first diagram indicate that the iterated OPEP was subtracted. 
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FIG. 2. Variable phases (normalized to 1) for the wave at the energies E = 40 MeV (dotted 
line), E = 100 MeV (dashed line) and E = 300 MeV (continuous line). 
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FIG. 3. Variable phases for the waves 3 D% (dotted line), 3 Gs (dashed line) and the mixing 
parameter €3 (solid line) for background potentials cut at 0.7 fm (left) and 1.0 fm (right). 
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FIG. 4. Phase shifts for the channel (T = 1,5 = 0) due to the OPEP cut at 1.0 fm (dotted 
line), the sum of the OPEP and the pure pion-nucleon TPEP cut at 1.0 fm (dashed line) and 
the full chiral potential cut at 1.5 fm and 2.5 fm (continuous lines labeled x(l-5) and x(2-5)). For 
comparative purposes, the predictions of the argonne vl4 potential are also included (dashed-dotted 
line). The experimental points represent the pn solution of SAID [23|. 



13 





-1 

-2 
-3 
-4 


-1 
-2 
-3 

-4 

0.0 
-0.5 
-1.0 
-1.5 




" x(2-5) 

, I , , , , I , , , ,_J , L 

100 200 300 




100 200 300 

Energy (MeV) 



w 

CD 
CD 

CD 



w 

CD 

W 

cd 



2.5 
2.0 
1.5 
1.0 
0.5 

0.0 

-0.5 

-1.0 

-1.5 

-2.0 
0.8 

0.6 

0.4 

0.2 

0.0 



, , , 1 , 

M 


, | , , 
(c) 


i 








i 




- 50 


100 


150 


200- 


, , , i 






n , . 

•^(1.5) - 
' " x(2-5) 


100 


200 




300 - 


; 3 h 4 






x(l-5) 

















l , 



3 100 200 300 

Energy (MeV) 




FIG. 5. Phase shifts for the channel (T = 1, S = 1); (a) uncoupled waves; (b) J = 2 coupled 
channel; (c) J = 4 coupled channel; (d) J = 6 coupled channel. Conventions are given in Fig. 4. 
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FIG. 6. Phase shifts for the channel (T = 0, S = 0); conventions are given in Fig. 4. 
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FIG. 7. Phase shifts for the channel (T = 1, S = 1); (a) uncoupled waves; (b) J = 3 coupled 
channel; (c) J = 5 coupled channel. Conventions are given in Fig. 4. 
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